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Abstract 

The production of J/^> 's in nuclei is re-examined and the critical role of Feynman 
x distributions in the study of photoproduction and hadronic production of both J/^'s 
and Drell-Yan dimuons is expatiated. The need to consider both initial and final state 
interactions and the key effect of initial state interactions on the total J/^f cross section 
are demonstrated. 

We first present a theoretical study of the expected functional form for the pro- 
jectile and target atomic number (A,B) dependence of the J/ty suppression, S, and 
demonstrate why it cannot fall exactly exponentially with A 1 /' 3 + B 1 ^ 3 , since such a 
term is multiplied by a weak enhancement factor. We use a Woods-Saxon Monte Carlo 
simulation to obtain distributions of the numbers of collisions prior and subsequent to 
the charmonium production. With those results we relate the mean number of nucleon- 
nucleon collisions, < n > to A 1 / 3 + B 1 / 3 , and give a new analytic functional form for 
the A dependence of the suppression. 

Finally we carry out a full Monte Carlo calculation of S, including an initial state 
prior energy loss in keeping with the measured Feynman xp distributions in tc-A and 
p-A interactions. We use an energy loss parameterization which is consistent with the 
energy loss from minimum bias data fitted in ISAJET. We also use an open charm 
absorption cross section, a oc , constrained from photoproduction of J/^'s , to analyze 
the latest data on the J/\I/ suppression. 

The role of color screening, the time evolution of off-shell charmonium to the on-shell 
J/\I/ , and the energy dependence of the charmonium-nucleon inelastic and absorption 
cross sections are discussed. A program of related measurements, necessary for better 
understanding the J/^f suppression, is presented. 

We find no anomalous suppression when comparing our results to published data. 



1 



I. Review 



Shortly after the proposal U of the possibility of seeing evidence for deconfinement 
of charmonium in nuclear reactions, use was made of studies of soft processes in rela- 
tivistic nuclear reactions to examine the rich Feynman x distributions of J/ty 's that 
had been measured in production of J/\I/'s in pion and proton interactions in nuclei. 
From the outset [[J, charmonium absorption, as the only source of J/ty suppression, 
was ruled out by a glance at the many examples of the strong xp dependence of the 
p-A/p-p and 7r-A/7r-p ratios of the J/\& yields. In that first work, only the energy 
loss of incoming hadrons before the charmonium production and energy loss due to 
rescattering of charmonium off nucleons in the final state were considered. A more 
careful examination of many examples of xp distributions for both Jfty and Drell-Yan 
production was then carried out [Q. 

A direct determination of the charmonium absorption cross section, a oc , from pho- 
toproduction of J/^f 's was then obtained ||. In photoproduction there are no initial 
state interaction effects to confuse the interpretation, as in production from hadron 
probes. We obtained the result, a oc = 6.6 +/- 2.2 mb. 

It is important to understand that there is, presently, no way to calculate the 
non-perturbative effect on the nucleon structure functions resulting from prior soft 
scatters before the high Q 2 production of vector mesons. Thus the "energy loss" must 
be treated as a parameter that can be obtained, in principle, from analysis of the 
Feynman x distributions. 

To carry out these early calculations we used the well-known ISAJET model of 
Frank Paige to arrive at an expression for the approximate functional form of the i/i 
dependence of the energy loss in soft collisions prior to the J/\I/ production. The energy 
loss per collision is found in low pt interactions to vary approximately as dyfs/dn = a 
+b-y/s . The functional form is useful but the values of a and b obtained from ISAJET 
are only a guide. For many calculations of the xp distributions it is sufficient to use 
dy/s/dn = a alone. We studied d^/s/dn in the range .4 to 1.0 in our early work. (Ref. 

i 

In Ref. || , we presented two key examples of our earlier calculations of the Feynman 
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x distributions for a pion- and a proton-induced reaction. Two samples of data are 
reproduced in Fig. 1, since they clearly illustrate the many effects that are responsible 
for J/^f suppression: 



150 GeV 71+ Pt/p -*J/\|/ 
1 .0 i ' ' ' ' i ' ' i ' i ' ' i ' i 




0.0 0.2 0.4 0.6 0.8 1.0 



800 GeV p + W/D-»- J/y 
1.0 i 




x F 



93-0028a,b (comb)new 

Figure 1: 

1) Loss of J/^> 's due to absorption: It is generally assumed, since there are no data 
to test the hypothesis, that a oc is energy independent. Thus the charmonium absorption 
alone can only reduce S(xf) independent of xf ■ A flat reduction independent of xf is 
not what the data show. All measured reactions show falloff of the ratio of the nuclear 
to nucleon production with xp ■ Further, suppressions at x F = are not the same for 
pion and proton induced reactions, additional proof that final state absorption cannot 
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be the only mechanism for the suppression. 

2) Initial state energy loss: This has two effects: a) the yield is reduced because 
of the exponential dependence of the J/^f yields on 1/y/s , e~ M /^, and b) events are 
shifted to lower xf ■ Both effects appear in all the measured xf distributions. 

3) The outgoing charmonium interacts with nucleons inelastically, losing further 
energy by soft collisions: J/^ + n^J/^ + X. This will affect the shape of the x F 
distributions as well, forcing events to lower xp ■ 

(In Fig. 1 the absorption cross section was taken as 9 mb and d^fs/dn taken as 
.4 but a smaller cross section and somewhat larger dy/s/dn would provide a similar 
fit. Also, an energy loss of the outgoing J/^> of AE/E = .1 resulted in a need for an 
inelastic cross section, a ine i, of about 11 mb.) 

Clearly the Feynman x distributions are complicated, but they are rich in informa- 
tion. Fortunately the inelastic J/\I/ interactions will not affect the total J/ty suppres- 
sion ratios. (However, they will affect the rapidity distributions and thus affect the 
rapidity acceptance, producing a possible experimental A dependent effect on S.) 
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II. Bases for Suppression Calculations 



We digress to emphasize some of the constraints that time dilation places on mul- 
tiple scatterings in nuclei: Immediately after a low p t proton collision on a nucleon in 
a nucleus, the proton is in an "excited" state. It is " off-shell" . The asymptotic final 
state of the proton is a real "on-shell" proton. Real pions and other produced mesons 
appear later in the asymptotic state. The time for this evolution in the proton center 
of mass is usually assumed to be roughly set by the uncertainty relation, AEAt = H. 
t is then time dilated so that the final "asymptotic" state of on-shell proton and pions 
occurs much later. 

(One would like to believe this argument is quantitative, but the time evolution is 
not presently understood, as recent studies of color screening have indicated || . 

For light nuclei it is known that the hadronization occurs when the hadrons have 
left the nucleus 0. Since the measured multiplicity distributions in p-D, p-He, and He- 
He interactions are well fit from the p-p multiplicity distributions and the probabilities 
of making n wounded nucleons, neglecting any final state effects is reasonable for light 
nuclei. However there may be final state corrections needed for the multiplicity or E t 
distributions in heavy nuclei. Thus, one cannot use the p-p measured E t distributions 
directly in predicting the nuclear E t distributions. 

In this picture it is the excited off-shell proton, p*, which collides with nucleons 
to produce the charmonium in the following scatterings. We have made, in all our 
calculations, the assumption that the p* - nucleon cross section is little different than 
the p-nucleon cross section, on the possible ground that high energy nucleon-nucleon 
cross-sections are essentially geometric, as the energy independence of the p-p cross 
sections suggest. However, if the cross sections were to increase as the protons become 
excited, calculations of the number of soft collisions, prior and subsequent to the char- 
monium production, would be enhanced, increasing the energy loss effect. It is useful 
to appreciate this difference in the low p t collisions in p-A and A-B collisions. 

However, in A-B interactions, the charmonium may be interacting with nucleons 
struck by other incoming nucleons, so the A-B absorption is not on ground state nucle- 
ons but sometimes on excited nucleons. This too is ignored in the model but it points 
out that J/\I/ production in p-A reactions cannot strictly be treated identically with 



5 



A-B interactions. 

The high Q 2 charmonium state, however, has a short lifetime, so that At is much 
smaller than in low pt collisions. Therefore it is believed to be mainly made within the 
nucleus where it can interact with nucleons further along the path and be absorbed by 
reactions leading to open charmed particles. 

Because the pions and kaons made in the soft interactions mainly materialize outside 
the nuclei, these "co-movers" should not appreciably affect the suppression. We do not 
consider them and we shall see that they are not needed to account for the measured 
suppressions, although they might have a small effect on S. (Aspects of this subject 
are treated theoretically, and data on pion production by muons has been examined 
in a study of the x dependence of pion production in nuclei.) "A Dependence of 
Hadron production in Inelastic Muon Scattering and Dimuon Production by Protons" , 
S. Frankel and W. Frati, Phys. Rev. D51 4783 (1995). 

It is clear that the soft energy loss of the pertinent constituents, gluons and quarks, 
which produce J/^> 's and Drell-Yan pairs, and which should be different for proton and 
pion projectiles, will have to be parameters obtained from the Feynman x distributions. 

In our work Ref. H, studying disjoin for J/ty 's, we observed empirically that the 
value of d^/s/dn was smaller, = .2, for Drell-Yan pairs. This latter value is very crude 
but provided a rough fit to our analyses of various early experiments 0. 

As shown in Fig 1, reasonable fits to both 800 GeV proton data and the 150 Gev 
pion data with a charmonium- nucleon inelastic cross section of about 11 millibarns 
were obtained, assuming no color transparency effects were present. (It was assumed 
in these analyses that no color screening was present and that o oc was xf independent.) 

It was abundantly clear that energy loss in the initial state as well as the absorption 
into open charm played important roles. The initial state energy loss was crucial to 
the J/ty analyses and was small but not zero in the Drell-Yan data. 

Thus, in this paper we return to examining the new data on J/ty suppression in 
heavy nuclei since the recent analyses [TO] make assumptions which the Feynman x 
data contradict. 
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III. Analytic Study of the A Dependence of the J/ty Suppression 



In this section we see what insights can be obtained from collision of a 2 lines of 
interacting nucleons to obtain simple formulas for the effects of energy loss and open 
charm absorption on J/\I/ and dimuon production. As shown in Fig. 2. N and M 
are the number of nucleons in each nucleus along a collision line and n and m are the 
nucleon positions in the line. In a nucleus-nucleus collision there will be a distribution 
of such line-on-line collisions and the final result will require averaging over the N + M 
distribution. We shall describe such Monte Carlo results later in this paper but the 
analytic calculation for the line-on-line expression for the J/ty suppression contains 
interesting information. 

N=6 M=7 



— ► oooooo ooooooo -<-- 

654321 1234567 

t t 

97-0231 n m 

Figure 2: 

We shall assume for this calculation that all the measurements are taken at the 
same projectile energy so that e~ a represents the reduction in yield per soft collision 
due to energy loss prior to the J/\I/ production and e _/3 represents the absorption loss 
in a single J/\I/ -nucleon collision subsequent to the production. (No difference in any 
of the cross sections due to proton-neutron differences appears in this simple analysis, 
although it will affect yields at high A 1/s + B 1 / 3 .) 

Consider the J/^f production in the scattering of nucleon n = 2 on nucleon m = 3 
as shown in Fig. 2: n = 2 is slowed down by soft interactions with m — 1 and 
m — 2, while m = 3 has been slowed by interaction with n — 1. Our calculation sums 
over the n, m variables. Similarly, nucleons 3-6 in iV and 4-7 in M interact with the 
produced charmonium, these collisions absorbing the charmonium into free charmed 
particles, with an absorption cross-section denoted as a oc . Because of the high Q 2 
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(low cross-section) nature of charmonium production, the number of J/^'s produced 
is proportional to NM, all nucleons being equally likely to interact to make a J/f . 

The energy loss effect comes from the strong \fs dependence of the J/\l/ cross section 
which is known to be given by: 

1) <7j/« = e-< M l^ 

with 7 = 14.5 and M = 3.1 GeV. 

It is important to note that this expression is yfs dependent and the measurements 
cover a wide variety of \fs . 

Taking into account prior energy loss, this becomes: 

2) e -jM/(^o- An ) = e -lM/(^ x[l-An/^ ]) 

where A = d^fsjdn and n is the number of collisions prior to the J/ty production.. 
Expanding the denominator, it can be approximated by e -^ M l^)( 1+An l s °) . 
Thus we have the simple linear approximation: 

2') e — an = e~ ('yMA/s )n 

(Note that the energy loss effect depends on s and not y/s as in the cross sec- 
tion for J/\I/ production and that the lower the energy of the reaction, the larger the 
suppression.) 

As we shall see, this is a fair approximation and we shall use it in our analytic stud- 
ies. But we shall use the exact formulation in our Monte Carlo calculations reported 
later in this paper. 

Thus, for the interacting nucleons shown in Fig. 2, we have: 

3) energy loss: e -^[(n-i)+(m-i)] 

Similarly we can parameterize the loss into open charm where j3 is proportional to 
the open charm cross section, aoc, assuming that ooc is energy independent. 

4) open charm: e -Pl(N-n)+(M-m)} 
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This product can be re-factored as 

open charm: = e -P[{N-i)HM-i)} x e +/3[(n-i)+(m-i)] 



5) Including both energy loss and open charm we obtain: 

Sp = e -/3[(N-l)+(M-l)] x e -(a- /3 )[(n-l)+(m-l)] x 

Summing over the nucleons, this can be rewritten as: 

6) S p = e -^ +M ) x Ef Ef e -(«-W*-i) e -(°-0)CJ-i) 

Carrying out the sums and rearranging, we get the final result: 

7^ q _ „-(cH-/3)(iV-l+M-l)/2 v \ sinh(a-f3)N/2 ir sinh(a-(3)M/2 i 
'J°p- e x \-Nsinh(a-p)/2i\-Msinh(a-j3)/2i 

(N-l + M-l)/2 can be considered the mean number of prior or subsequent collisions 
in the row. We shall see later how this is related to the mean number of collisions 
calculated from a Woods-Saxon distribution which we relate to A 1 ! 3 + B 1 ^ 3 . 

We denote the bracketed factor in 7) as C . Since it increases with iV and M, C is 
an enhancement factor as opposed to the exponential suppression factor it multiplies. 

Note that the cancellation of a and (3 in C helps in making C small and the 
exponential fall-off a better approximation. 

We will return to this effect later when we discuss central A — B collisions, where the 
enhancement effect can be seen to be an appreciable contribution to the A dependence 
of S. 

A useful approximation for C is: 

8) C ^ [1 + (a - (3) 2 (N 2 - 1)/6][1 + (a - (3f{M 2 - l)/6] 

Thus we have demonstrated that the cross section does indeed fall off exponentially 
with a + P except for the enhancement factor, C, given by the sinh expression. The 
enhancement effect is understood by recognizing that the interactions at the very front 
of the nuclei have only open charm absorption while the interactions at the very rear 
of the nuclei have only energy loss. However, at other positions the effects are both 
present. 
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We now turn to a comparison of J production to Drell-Yan production. In this 
case there is no open charm contribution for the Drell-Yan production, only energy loss 
entering into the ratio, denoted as Sdy The result is: 



9) S dy = e -(a-a')(iV-l+M-l)/2 x Q ^ where 

r t sinh(a-P)N/2 -\T sinh(a-P)M/2 -\ / r sinh(a')N/2 i r sinh(a'M/2 

^dy- iNsinh(a-P)/2i\-Msinh{a-P)/2i / lNsinh(a'/2)l\-Msinh(a'/2) 



In this expression a' is the Drell-Yan energy loss parameter which is considerably 
smaller than that for J/ty production. We have estimated dsfsjdn = 0.2 in our prior 
studies of Drell-Yan production in nuclei. |4j [[| 

Prediction: The "slopes" of the logarithmic plots of S p vs the mean number of 
scatters should be somewhat larger than that for Sdy 

We now turn to an extension of the analytic calculation, which depends on the mean 
number of scatters in a line of interactions, to integration of the line interactions over 
the spatial configurations of the nuclei with the goal of relating (N — l)/2 + (M — l)/2 
to A 1 / 3 + I? 1 / 3 . We have shown that the J/^f suppression that would take place if 
a string of N nucleons interacted with a string of M nucleons along a line is given 
approximately by 

10) S = e-(°*w * 1 

The C of equations 7, 8 is a function of a , (3 , N and M but is close to unity so 
we will ignore it for the present discussion. In this form f3 represents the loss due to 
J/\l/ absorption but a represents an approximation of the effect of prior energy loss of 
the incoming nucleons on the J/^f yield. 

The question we address here is how ( JV ~ 1 )+( M ~ 1 ) depends on the quantity K(A 1 ^ 3 + 
B 1 ! 3 — 2) since this is the variable that many workers use in plotting the suppression 
data. (We choose to keep the -2 in these expressions since, for p-p collisions, these 
quantities are zero, so S is unity.) 

It is clear that equation 10) would have to be integrated over the probability that 
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one gets N and M nucleons in a line, averaged over the impact parameter in a nucleus- 
nucleus collision. 



Another question is whether the needed average is some simple function of (A 1 / 3 + 
5V3_2). 

To answer these questions we do not (at first) go to a full scale Monte Carlo cal- 
culation of the J/\I/ suppression but first consider only the distribution of prior and 
subsequent collisions, n prior = n su b s = n, before and after the J/\l/ production. 

These are needed to examine eq. 7. We obtained these distributions from a full 
Monte Carlo, using Woods-Saxon representation of the nuclei. 

The averages for p-A and A- A interactions are shown in Fig. 3. In this figure 
n-A + nB is the sum of the prior scatters in both nuclei. 




Figure 3: 

(For deuterium and helium we have used the best spatial distributions, as was done 
in our search for deconfinement in p, d, and a experiments at the ISR citeakesson. 
This shows up in the slight curvatures at very low A.) 
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The calculations show that the slope of < n > vs (A 1 / 3 + B 1 / 3 ) is essentially the 
same for p-A and A-A interactions, and a little reflection will allow the reader to realize 
that this is what is to be expected. We can then ask what value of K is needed to 
calculate < n >. What we find from the slopes in the figure is that K = .46. 

The fits to a straight line, < n > = a +b(A 1 ^ 3 + B 1 ! 3 -2), are given in Table 1. 

Table 1. 



b 



p-A average (prior) .46 .05 

A-A average (prior) .46 .00 

p-A central (prior) .63 .00 

A-A central (prior) .51 -.10 



We also note from the curves for central, i. e., zero impact parameter, collisions, 
Fig 4, that for a central AB collision the values of < n > are somewhat different than 
those of an AA collision with the same (A 1 / 3 + B 1 ^ 3 ) value. However, as expected, 
A — B and A — A average collisions should fall on the A — A curve at the appropriate 
A 1/3 + B 1 ' 3 . 

One might guess that a good approximation for the averaging over all n might be 
to replace (N — l)/2 by < n > in eq. 7. Instead, we will carry out the averaging to 
test this hypothesis. 

If one examines the distribution of prior collisions, n, when the average number is 
< n >, the normalized distribution is given approximately 0] by: 



11) F(n) = e~ n / <n> / <n> 

Therefore we can average S(n) over the n distribution to obtain S(< n >) 
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Figure 4: 

12) S{< n >) = j e <*n e -n/<n>i <n> dn 
This is just 

13) S{< n >) = ( a+1 / < ^ > ) <n> = i+a <n > 
Using 

14) 1/S(< n>) = l + a<n> 

To verify this formula, we can plot l/S(< n >) vs a, which should then be a 
straight line whose slope is the mean number of collisions. To do this we have used our 
complete Monte Carlo calculation for S(ct) for average Pb-Pb collisions to determine 
the validity of the approximate calculation of eq. 14. We have chosen a's corresponding 
with hypothetical absorption cross-sections from 4.08 to 9.4 mb. 

Fig 5 shows this plot. It is indeed linear with a slope of 0.44. which is, within error, 
the same as the value determined from the Monte Carlo calculation of < n > for a 
Pb-Pb collision, shown in Fig 3. 

Thus the theoretical assumption for an exponential shape of the n distribution is 
apparently a good approximation. 



13 




Figure 5: 



As a second check, we can examine, for fixed a, how the Monte Carlo calculation of 
< n > and the approximate calculation of < n > from eq. 14 depend on A 1 ^ 3 + B 1 ^ 3 — 2 
for p-A collisions. This plot is shown in Fig. 6 showing that both methods agree for a 
range of A's from P-C to p-U. 

It is very useful, therefore, to note that one can use the analytical expression of 
equation 1) and the parameter K = .46 to calculate effects of absorption whether they 
are due to disappearance of the J/\l/ or due to an exponential decrease in yield as the 
result of prior energy loss. 

Thus we conclude that an effective L for plotting the suppression should be 

15) L eff = AQ{A 1 ' 3 + B 1 / 3 - 2) 
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Figure 6: 



IV. Complete Monte Carlo Calculations of the J/ty Suppression 
A. Charmonium Absorption Only 

We first turn to a calculation of the suppression on the assumption that initial 
state energy loss can be completely neglected. Fig 7 shows the effect of varying the 
absorption cross section obtained from our full Monte Carlo and based on the same 



Woods-Saxon distribution used by M Nardi [13|. 

/,From this figure we see that the cross section that fits the data is greater than 9 
mb and thus is in disagreement with the cross section obtained from photoproduction. 
This alone points to the need for energy loss. However Ref. |IB| has done what is 
presumably the same calculation and obtain a cross section of 7.6 mb. We believe 
that this might be due to the following difference in the details of the calculation: In 
calculating absorption effects it is important to recognize that the nucleus is composed 
of finite nucleons and that there is a counting problem when one attempts calculations 
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Figure 7: 



using a smooth density distribution. If one integrates along a path up to a point within 
the nucleus and then asks for the probability of a collision on the way out of the nucleus, 
the struck nucleon is included in the absorption path, overestimating the absorption 
and thus requiring a smaller absorption cross section. One cannot have the struck 
nucleon absorb the charmonium without including this effect in the p-p cross section 
as well. This effect also occurs in the calculation of the effects of color screening in high 
p t e-p and p-p collisions. There the overcounting resulted in a claim for observation of 
color screening. The difference in the two methods is easily demonstrated ||. 
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B. Energy Loss Only 



1.0 
0.9 

0.8 

_ 0.7 

CL 

B 0.6 
D 

X 0.5 
CO 

I 

< 

ll 

DC 

0.3 









Energ 


y Loss Onl 


/ 






j*>^^r^^^^ 















j«? 




/- — 




38.8 
29.0 









□ 






17.4 




o 


Vs =38.8 


♦ 


O 

□ 








u- 

♦ 


•/s -29.0 

Vs =17.4 

a = 6.3 mb 

oc ; 




□ 

♦ 

♦ 


O 






d 

dn 


= 0.5 +0.018 Vs 






C 







i i i 


1 1 1 il 1 1 1 


1 1 1 


1 1 1 


♦ 

1 1 1 


i i i 





"0 2 4 6 8 10 12 14 

A 1/3 + B 1/3 



Figure 8: 



In this subsection we show the effects of energy loss with and without absorption. 
Fig 8 shows a plot of average collisions for several values of yfs , corresponding with 
ones used in present experiments. The upper lines show the effects of energy loss alone. 
The lower points show the effects of changing y/s when both absorption and energy 
loss are computed. 
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C. Combining Absorption and Energy Loss 
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Figure 9: 



In this final section we turn to a complete study of both the energy loss and absorp- 
tion effects without the exponential energy loss approximation. It is obvious from eqs 
1) and 2) that energy loss will produce deviations from pure exponential behavior of S 
vs A 1 / 3 + B 1 / 3 . If the energy loss were given by the approximate relation, namely e~ an , 

we have seen that S falls off almost exponentially with a. However, as we have seen 
previously, the actual energy dependence is more complicated. Unlike the so-called 
"scaling factor" of eq. 1, which depends on y/s , the energy loss factor is more sensitive 
and depends on s, as shown in eq. 2. The effect of using the true dependence of the 
J/\P cross section on energy is to introduce curvature into the plot of S vs the mean 
number of collisions. The direction of the curvature is to depress S at high A and low 
y/s which, as we shall see from our Monte Carlo calculations is in agreement with what 
the data show. 
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Fig 9 shows the results using equation 2) with the parameters d^/s/dn = .5 + .018 
and a oc = 6.3 and 7.9 mb. These are reasonable choices of the parameters and 
are consistent with the photoproduction data on o oc . However, other choices, within 
limits, give equivalent results. It is worth pointing out, however, that the data are taken 
at different energies and with different experimental arrangements and no systematic 
errors have been provided. We have marked the different points so the reader can 
decide whether slightly larger a oc with a smaller dyfs/dn will give a "better fit", but 
we will not do so, in view of the need for three parameters which are not accurately 
known. 

This simple counting model, incorporating energy loss, accounts for the present 
data without the need to invoke any "new physics" . 
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VI. Central Collisions 
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Figure 10: 



We have seen in section II. that there is an enhancement factor C that multiplies 
the exponential dependence S on A 1 ! 3 + B 1 ^ 3 . This was shown for the special case of 
interactions of two lines of nucleons. One can verify this, of course, using the full exact 
Monte Carlo. The largest number of scatters would occur in measurements that trigger 
on high transverse energy, so this type of measurement would find deviations from true 
exponential behavior. Since high E t triggers favor central, zero impact parameter, 
collisions, we choose to examine central A-A collisions. Fig 10 shows a plot of central 
A — A collisions for a oc = 7.9 mb. The energy loss is set equal to zero to demonstrate 
the effect. It is clear from the full MC calculation that S turns up at large A. Once 
again we see how simple exponential extrapolations to large A can be misleading. 
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VII. Examination of the y/s Approximation 
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In section III. we carried out an analytic calculation based on the approximation in 
the energy loss equation, going from eq 2 to eq 2'. To see this effect in more detail, we 
show in Fig. 11a comparison of the two results. It appears that the approximation is 
noticeable but produces only a small effect on the suppression even at high A 1/3 + B 1 ^ 
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VIII. Discussion and Future Analyses 



Because the soft processes that take place in nuclei are not able to be handled by 
perturbative QCD, it is essential to examine how J 's are made in a large variety of 
experiments. It is especially important to study the energy loss effect since it appears 
in all particle production in nuclei. 

1) Photon and electron probes, having no initial state interactions, should be used 
to determine both the xp dependence of the open charm cross section and the energy 
loss of the outgoing J/ty . Both are presently unknown, as is the J/^f -nucleon inelastic 
cross section. Several incident energies will be needed to unravel these three quantities. 

2) Feynman x distributions of Drell-Yan dimuons are poorly known. This reaction 
studies only the initial state interactions and is the companion to the photo-production 
experiments. The Drell-Yan data are the only data that determine d\fs/dn uniquely, 
since there are no final state interactions. This information is needed in comparing J 
production relative to Drell-Yan production. The Feynman x distributions at xf = 1 
are constrained by a knowledge of the Glauber coefficients; in fact, R is given by the 
probability that the pair is made in the first collision. Thus the xf = 1 measurements 
are a check on Glauber coefficients and the accuracy of the nucleon spatial distributions. 
The Xp A dependence directly checks this information, which enters into every nuclear 
reaction. 

3) Finally, the Feynman x distributions for the nuclear J/ty cross section relative 
to p-p and relative to the Drell-Yan can be examined and compared for internal con- 
sistency. 

Even these measurements may not allow one to get a full understanding of the 
processes for the following reasons: 

1) As pointed out by Hufner et al. |14]], that part of charmonium that is in the 



singlet state should be color screened and therefore a oc can easily be xp dependent. 
It cannot be a major effect, since the fall in the ratio K(xp ) with xp seems to be 
dominated by energy loss, but any color screening will contribute to an increase in the 
values of S. 

2) The time evolution of the off-shell charmonium to the on-shell J/^> configuration 
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is poorly known. There are, in fact, no estimates other than AEAt arguments and no 
calculations of the time evolution of the transition of a c—c color-screened configuration 
to the J /ty on-shell configuration of charmed quarks and gluons, especially since this 
is surely dependent on the density of gluons in normal nuclear matter. This must be 
known since the final state interactions of charmonium and the J /iff are not necessarily 
identical. 

We conclude that careful study of J /iff production in nuclei can teach us about the 
interplay between the non-perturbative and perturbative processes of the J /iff , but 
that claims for "new physics" from total yields have little bearing on these interesting 
problems that, so far, are poorly understood. 

IX. Conclusions 

1) Analytic examination of the J /iff suppression, S, shows that S cannot the be 
exactly exponential in A 1 / 3 + B 1 ^ 3 and therefore that such an extrapolation to large 
A 1/3 + B 1/3 is invalid. 

2) All Feynman xp p-A data demonstrate that final state absorption of charmonium 
cannot be the only contributor to the suppression, so calculations of suppression in A-B 
interactions must include the effects of energy loss in soft scatters prior to charmonium 
production. 

3) Complete calculations using a Woods-Saxon nucleon distribution shows that 
there appears to be nothing "anomalous" in the suppression observed in Pb-Pb colli- 
sions. 
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Figures 



Fig. 1 Measured Feynman x distributions, for the nucleus to proton target ratio, 
R, of J/\I/ production. (Also shown are fits to the data, including a crude estimate of 
the J/ty - nucleon inelastic cross sectiom in the final state taken from Reference ||. 

Fig. 2 Schematic of interactions of nucleons in a line: representation of the collision 
of a line of N nucleons with a line of M nucleons. 

Fig. 3 Results of Monte Carlo Woods-Saxon calculation of the mean number of 
collisions prior to (or succeeding) an interaction producing charmonium in p-A and 
A-A average collisions. 

Fig. 4 Results of Monte Carlo Woods-Saxon calculation of the mean number of 
collisions prior to (or succeeding) an interaction producing charmonium in A— A central 
collisions of zero impact parameter as well as for average A — A collisions. Several A — B 
collisions show departures from the straight line slopes, but only for central collisions. 

Fig. 5 Plot of 1/S = a+ < n > from the exact full Monte-Carlo calculation, using 
the values of < n > for Pb-Pb and varying a oc from 4.1 to 9.4 mb. This plot is used 
to demonstrate the validity of this new equation which has been derived using the n 
distribution of equation 11. 

Fig. 6 Comparison of the A 1 ^ 3 + B 1 / 3 dependence of the mean number of scatters 

< n > m c, taken from the full exact Monte Carlo calculation of S with the value, 

< n > approx, using the derived formula, 1/S = a+ < n > approx, for the case o oc = 7.4 
mb. 

Fig. 7 Complete Monte Carlo calculation for pure absorption. The best fit disagrees 
with the value of a oc from photoproduction. 

Fig. 8 Final predicted values of the suppression for different values of a/s for a 
fixed a oc = 6.3 mb, showing how energy loss has different suppressions in experiments 
at different ^fs . 

Fig. 9 Comparison of Calculations with Data at Various Values of . The 
solid line shows the extrapolation of the lower yfs data. Calculations are shown for 
two values of o oc consistent with the photoproduction data for a reasonable value of 
dy/s/dn (See text.) 

Fig. 10 Calculated Suppression in Central A-A Collisions for a oc = 7.9 mb and no 
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energy loss, showing that the enhancement factor, C, shows up in the full Monte Carlo 
Calculation. 

Fig. 11 Effect of using the exact expression for the \fs dependence of the J/\I/ 
cross section, eq. 2 vs the use of eq 2'. The plot is shown for d^fsjdn = .4 + .012s and 
a oc = 6.9 mb. 



25 



References 



[1] T. Matsui and H. Satz, Phys. Lett. 178B 416 (1986) 

[2] "Predicting Ultrarelativistic Nuclear Interactions from p-p Data", S. Frankel and 
W. Frati, Nuclear Physics B308 (1988) 699 

[3] " J/\I/ and Dimuon Production in Ultrarelativistic Nuclear Collisions", S. Frankel 
and W. Frati, Quark Matter in Collision, World Scientific. (1988) 

[4] Dimuon and Vector Meson Production in Nuclei", S. Frankel and W. Frati , UPR 
50 pp unpub. UPR 0342 T (1990) 

[5] "The J/\I/ Absorption Cross Section in Nuclei and Color Screening", S. Frankel 
and W. Frati, Zeit . Phys. 57 225-7 (1993) 

[6] S. Frankel and W. Frati Physics letters, B291 368 (1992); S. Frankel, W. Frati, 
and N. R. Walet, Nucl. Phys A580 595 (1994) 

[7] Phys. Lett. 119B 464 (1982); Phys. Rev D28 2334 (1983); These were the ex- 
periments that verified, for multiplicity distributions, the validity of the Bialas, 
Bleszynski and Czyz hypothesis (Nucl. Phys. Bill 461 (1976) that the number 
of participants and not the number of scatters determined the event structure. It 
is now appreciated that the basis for the wounded nucleon counting is the time 
dilation effect. 

[8] "A Dependence of Hadron Production in Inelastic Muon Scattering and Dimuon 
Production by Protons", S. Frankel and W. Frati, Phys. Rev. D51 4783 (1995). 

[9] J. Badier et al. Zeit fur Phys. 20 201 (1983), S. Katsenevas et al. Phys. Rev. Lett. 
60 2121 (1987); A.S. Ito et al. Phys. Rev. D 23 604 (1981), P. Bordalo et al., Phys. 
Lett. 193 368 (1987), D. M. Aide et al., Phys. Rev. Lett. 64 2479 (1990). 



[10] See, for example, D. Kharzeev et al., |hep-ph/96122T7l 30 Nov 1996. 



[11] Akesson et al. Phys. Rev. Lett. 55 2535 (1985). 



26 



[12] Tabulation of Distributions of Scatters in Nucleus-Nucleus Collisions - UPR Re- 
port 776T, W. Frati and S. Frankel, 1997. 

[13] Private communication. We used the value of R, R = 1.19A 1 / 3 — 1.61A -1 / 3 , re- 
ported in D. Kharzeev et al. CERN-TH/96-328 BI-TP 96/53 (1996). 

[14] Phys. Lett. 258 465 (1991). 

jpfinal February 1, 2008 



27 



